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Popular Summary 


We derive routinely the concentration of aerosol particles, their properties and effect on 
solar radiation from satellite data, namely the MODIS instrument on the Terra platform. 
Aerosol is the haze part of urban pollution, smoke forest fires or desert dust. The aerosol 
information is derived over ocean and land. While ocean is uniformly dark and thus easy 
background to see the re (lection of sunlight from the aerosol, the problem is more 
complex over the land. MODIS detects sunlight reflected both by the land surface and by 
the aerosol. How to distinguish between them? 

We found that pollution and smoke aerosol is transparent at long invisible solar 
wavelengths. We therefore can use them to see how what the surface is composed of 
without aerosol interference, and try to use this information to guess how reflective the 
surface is in the visible part of the spectrum. Inferring the surface reflection of sunlight in 
the visible from the long wavelengths is based on collection of many measurements in 
different parts of the world that we did before the launch of MODIS. However we never 
had a clear physical understanding of this method. In this paper we for the first time use a 
model that simulates surface reflection across the solar spectrum for mixtures of soils and 
vegetation. We prove that our relationship between the surface properties in the visible 
and in the long wavelengths works most of the time, and found methods to correct them 
in cases it does not work. 
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Abstract: Data from the Moderate Resolution Imaging 
Spectroradiometer (MODIS) instrument that flies in polar 
orbit on the Terra platform, are used to derive the aerosol optical 
thickness and properties over land and ocean. The 
relationships between visible reflectance (at blue, 
and red, Pred ), and mid-infrared (at 2.1 /xm, p 2 ,) are used 
in the MODIS aerosol retrieval algorithm to derive global 
distribution of aerosols over the land. These relations have 
been established from a series of measurements indicating that 
Pbiue ~ 0.5 p red ~ 0.25 p 2 .i- Here wc use a model to describe the 
transfer of radiation through a vegetation canopy composed of 
randomly oriented leaves to assess the theoretical foundations 
for these relationships. Calculations for a wide range of leaf area 
indices and vegetation fractions si ow that p bluc is consistently 
about 1/4 of p 2 .i as used by MODIS for the whole range of 
analyzed cases, except for very darl soils, such as those found in 
bum scars. For its part, the ratio f r «/p 2 1 varies from less than 
the empirically derived value cl 1/2 for dense and dark 
vegetation, to more than 1/2 for bright mixture of soil and 
vegetation. This is in agreement with measurements over 
uniform dense vegetation, but not with measurements over 
mixed dark scenes. In the later case, the discrepancy is probably 
mitigated by shadows due to une< en canopy and terrain on a 
large scale. It is concluded that the value of this ratio should 
ideally be made dependent on ihe land cover type in the 
operational processing of MODIS data, especially over dense 
forests. 

1. Introduction 

A new generation of space-borne Earth observation 
instruments has been launched and can be used to monitor 
aerosols over land [King et ah, 1WJ. Satellite sensors differ in 
their measurement capability: MODIS on the Terra satellite is 
using a wide spectral range (0.47 to 2.1 /xm) to derive the 
aerosol optical thickness and properties [Kaufman etal., 1997a]. 
The Multi-angle Imaging SpectroRadiometer (MISR), also on 
Terra, retrieves aerosol properties using a modest spectral range 
(0.44 to 0.86 /xm) but observing the Earth from nine different 
directions simultaneously along the satellite track [Diner et al ., 
1989, 1998]. The POLarization and Directionality of the Earth’s 
Reflectances (POLDER) instrumLnt, operating in 1996-1997, 
using the same spectral range, observed the Earth from a wide 
range of angles in two dimensions and derived the aerosol 
properties with the help of polarization. The Along Track 
Scanning Radiometer (ATSR) instrument combines 4 spectral 
channels (0.55 to 1.65 /xm) with two view directions [ Veefkind, 
1998]. In all these instruments ttu measured signal consists of 
sunlight reflected by both the Eiirlh surface and the aerosol 
layer. Some assumptions about the surface reflectance properties 
and aerosol properties have therefore to be made to separate the 


aerosol contribution to the signal from that of the surface. Each 
instrument takes advantage of its unique configuration to derive 
useful information on atmospheric aerosols. 

In this paper, we concentrate on the MODIS technique 
based on the assumption that the ratios of the surface 
reflectances in the blue (0.47 /xm) and red (0.67 /xm) to the 
reflectance at 2.1 /xm are constant in space and time at the fixed 
values of 0.25 and 0.50, respectively. Our simulations were 
carried out at 0.44 /xm, which has similar characteristics to the 
MODIS 0.47 /xm channel. This ratio was established using 
remote sensing measurements in the United States [Kaufman et 
al, 1997b] and confirmed in an independent experiment in Israel 
[Kamieli et al, 2001]. The global application of the technique 
[Chu et at., 2002] shows that the method works over most 
terrains for which the surface reflectance at 2.1 /xm is smaller 
than 0.15. However no theoretical basis for this assumption had 
been developed so far. More specifically, the MODIS technique 
is based on the fact that sub-micrometer aerosol particles are 
transparent at 2.1 /xm (i.e., smoke and pollution but not dust). 
The MODIS algorithm first considers the reflectance at 2.1 /xm 
at the top of the atmosphere, and then uses this reflectance value 
to estimate the surface reflectance in the blue and red channels, 
using the empirical relationships described above. The aerosol 
optical thickness and properties are derived directly from the 
measured excess reflectance in the blue and red channels. In this 
paper, we provide some theoretical basis in support of the fixed 
spectral ratio approach, using surface bidirectional reflectance 
model of Gobron et al [1997]. 

2 Model simulations 

We simulate the angular and spectral properties of fully and 
partially vegetated surface using the plant canopy radiation 
transfer model developed by Gobron et al., [1997]. This model 
solves the radiation transfer problem in the case of finite-size 
oriented leaves uniformly distributed over a scattering soil. The 
leaves are not clumped into trees, but are distributed in the space 
as illustrated in Fig. 1. The model computes the Bidirectional 
Reflectance Factors (BRFs), p™^, for mixed surfaces (see fig. 
1) with fraction f, of soil and f v of vegetation (f $ +f v =l) using a 
linear mixing between them: 

Pmixcd (a 0 ,a) = S i f iPl (a 0 ,fi) (D 

Here, f, is the relative fraction of the pixel area occupied by 
surface type i. Allowed surface types include either vegetation 
with a given leaf area index (LAI) or soils, while and Q 
denote the solar illumination and satellite observation directions, 
respectively. The model simulations are limited here to one type 
of leaf properties and three types of soils. Table 1 summarizes 
the leaf and soil spectral properties. The latter were extracted 
from the database of Price [1995] and represent a wide range of 
actual soils. The leaf spectral properties are based on the leaf 



model of Jacquemoud and Baret [1990] for typical standard 
conditions of chlorophyll concentr ation and water content in a 
healthy leaf. It is important to note that the optical properties of 
the leaf and soils used in the present calculations were selected 
well in advance and independently of the MODIS algorithm and 
the measurements that serve the basis for the surface spectral 
properties in that algorithm. The radiative transfer model has 
been extensively evaluated and tested against a set of three- 
dimensional models [Pinty et at. , ' 001]. The linear assumption 
used in eq.l used to represent the angular reflectance of mixed 
pixels is assumed accurate enough for medium to low spatial 
resolution sensors. The calculations are performed for a wide 
range of leaf area index (0.05 to 5l, vegetation fraction (0.25 to 
1), view zenith angle (0° to 40°) ;ind solar zenith angle (1° to 
51°). The relative azimuth between the Sun and the observer is 
0°, 90° or 180°. 



Figure 1: Artist view of the stmcturally homogeneous leaf 
distribution implemented by the canopy radiation transfer model 
(top), and of the mixed surfaces used in our simulations, here for 
a fractional vegetation cover of 50% (bottom). 

Table 1: Spectral properties of leaves and soils used in the 
simulations. The leaf spectra are based on the leaf model of 
Jacquemoud and Baret [1990] for typical standard conditions of 
chlorophyll concentration and water content in a healthy leaf. 
The soil spectra are derived from the database of Price [1995] 
and represent a wide range of actual soils. 


Reflective surface 

wavelength 

/ 0 44 pm 

0.67 pm 

2.1 pm 

Leaf reflectance 

0)54 

0.054 

0.148 

Leaf transmittance 

0 X)5 

0.014 

0.172 

Dark soil reflectance 

0 068 

0.126 

0.140 

Medium soil reflectance 

0 102 

0.238 

0.357 

Bright soil reflectance 

0 139 

0.351 

0.536 


3. Results 

The plant canopy radiation transfer model [Gobron et a/., 
1997] is used to calculate the spectral reflectance of mixed 
scenes using eq. 1, where f v represents the fractional vegetation 
cover (f v = 0.25, 0.5, 0.75 and 1.0) and (1- fj is the fractional 
soil cover. The leaf area index (LAI) varies from 0.05 to 5. 
Throughout the ensemble of simulations, the vegetation density 
(f v *LAI) thus varies from 0.01 to 5. Note that shadows due to 
variability of the height of vegetation and terrain relief are not 
accounted for in these calculations. We shall evaluate the 
implications of this issue in the discussion section. Figure 2 
shows scatter plots of the simulated mixed surface reflectance at 
0.44 pm and 0.67 pm as a function of the reflectance at 2.1 pm 
for all cases corresponding to medium and bright soils. The 
vegetation density is represented in color, ranging from dark 
green for f v *LAI=5 to red for f v *LAI=0.01. Lines corresponding 
to the ratios of 1/4 and 1/2, for the blue and red respectively, as 
assumed in the MODIS algorithm [ Kaufman et al., 1997a] are 
superimposed. 
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Figure 2: Scatter plots of the simulated mixed surface 
reflectance at 0.44 pm (top) and 0.67 pm (bottom) as a function 
of the reflectance at 2.1 pm for medium and bright soils (see 
Table 1). The vegetation density is represented by the color bar, 
varying from dark green for f v *LAI=5 to red for f v *LAI=0.01 . f v 
is the vegetation fraction, and LAI is the leaf area index. Lines 
for the ratio of 1/4 and 1/2, assumed in the MODIS algorithm 
[ Kaufman et ai, 1997a] are also shown. The calculations are 
performed for view zenith angles (0° to 40°), solar zenith angle 
(1° to 51°) and azimuths of 0°, 90° or 180°. 



Quasi-linear relationships are observed between the blue 
(0.44 /un) or red (0.67 /cm) and the 2.1 /im channels for a wide 
range of simulated conditions. The simulated 0.44 to 2.1 / im 
ratio closely follows the assumed ratio of 1/4 over the whole 
range of surface reflectance values. In the case of the 0.67 over 
2.1 /im ratio, however, the simulated ratios are lower than 1/2 
for dark surfaces and larger than i/2 for bright surfaces. The 
figure shows a transition from ratios p re d /f P 2 .i around 0.25 for 
dense dark vegetation (typically forest), to ratios larger than 1/2 
for bright soils (see Table 1). In fact, measurements do show 
that for dense homogeneous forests the relationship between the 
reflectance in the 0.44, 0.67 and 2.1 /am is such that pbi uc ~ Prcd ~ 
0.25p 2 i [ Renter et ai, 2001; Gatebe et al., 2001], as predicted 
by the model. For mixed d irk scenes, however the 
measurements show that ~ 0.5 p rc d ~ 0-25 P 2.1 [fCaufman et 
ai, 1997a,b; Kamieli et al ., 2001]. 



Surface reflectance at 2,10 /j.m 

Figure 3: Similar calculations to those in Fig. 2, but including 
the dark soil reflectance (dots) from Table 1. The scale of the 
surface reflectance is limited here to a lower range of 0 to 0.2. 
The colors represent the leaf area index, LAI (not the f v *LAl 
value as in Fig. 2), which varies from 0.05 to 5. For each LAI 
value, the vegetation fraction varies from 0.25 to 1.0. The 
calculations are performed for view zenith angles (0° to 40 ), 
solar zenith angle (1° to 51°) and azimuths of 0°, 90° or 180 . 

Figure 3 shows similar calculations to those of Fig. 2, but 
this time including simulations with dark soils. The plots are for 
a limited reflectance range of 0 to 0.2 to emphasize the range 
that is used by the MODIS analysi s Note also that the color bar 
indicates the various LAI values instead of the vegetation 
density (f v *LAI). As already shown in Table 1, dark soils reduce 


the reflectance at 2.1 /mi much more than in the red or blue 
wavelengths, thus resulting in a much steeper slope. Such dark 
soil conditions are typical, for instance, of recent bum scars. 
Other soil conditions intermediate between dark and medium, 
would probably end up with a series of different slopes. 


It is interesting to see whether the ratios of Fig 2 depend on 
the view and solar illumination directions. The ratios of the blue 
and red reflectances to the 2.1 /m 1 reflectance are shown in Fig. 
4 for two solar zenith angles: 21° and 51°, representing 
conditions of high and low Sun. They are represented as a 
function of the view zenith angle. The figures correspond to full 
and half fractional vegetation covers, for LAI varying from 0.05 
to 5 and for conditions where the reflectance at 2.1 /m 1 is less 
than 0.25. 

100% Coverage 



07 

0.6 

0.5 

0.4 

0.3 

0.2 

■ 0.1 
0.7 
. 0.6 
0.5 
. 0.4 
. 0.3 
0.2 
■ 0.1 






-40 -20 0 20 40 

Vtswlng Angles (dsgr**s) 


-40 -20 0 20 40 

Visaing Angles (degrees) 


C u v m n" 


_ < 4ft 

i( -».ra 


* 

5 




5 4k ill 
5 *via 

*.10 

f fr.JKV 
^ *.30 
* ». 

•S^OA 



-*■: 0 sc 40 



~*J 0 50 *0 

V Wwi r i taflui 


Figure 4: Ratios of the blue (0.44 /im) and red (0.67 /mi) 
reflectances to the 2.1 /im reflectance, represented as a function 
of the view zenith angle for solar zenith angles of 21° and 51°, 
typical of high and low Sun. The figures are given for full (top) 
and half fractional vegetation cover (bottom), for conditions 
where the reflectance at 2.1 /im is less than 0.25. The LAI varies 
from 0.05 to 5 following the color convention used in Fig 3. The 
calculations are performed for view zenith angles (0° to 40°) and 
azimuths of 0°, 90° or 180°. 





It can be seen that the ratios ^e weakly dependent on the 
view zenith angle for low solar zenith angle. They are also 
weakly dependent on view zenith angle for vegetation mixed 
with soil- The dependence is stronger for higher solar zenith 
angles over full vegetation cover. MODIS will encounter high 
solar zenith angles at high latitudes e.g., over the boreal forests. 
In these cases MODIS analysis will underestimate the surface 
reflectance by 0.01 at 0.67 ym (negligible error at 0.47 ym). For 
these low solar zenith angles - 51 \ view angles -40° and full 
vegetation cover the corresponding underestimate in the aerosol 
optical thickness, At is relatively small. Using radiative transfer 
calculations for typical aerosol model, we get: At ~ 0.07. This 
error can be corrected by accounting, in the MODIS algorithm, 
for the angular dependence suggested by the model. 

4. Discussion and conclusions 

A model describing the transfer of radiation through a 
vegetation canopy composed of randomly oriented leaves is 
used to formulate a theoretical basis for the empirical 
relationship between the surface rellectances in the blue, red and 
2.1 gm channels assumed in the MODIS aerosol algorithm. The 
calculations were carried out for a wide range of leaf area index 
and vegetation fractions, using a linear mixing between 
vegetation and soil reflectances. 

We show that the reflectance n the blue is consistently 1/4 
of the reflectance at 2.1 ym for all analyzed cases over medium 
and bright soils. This finding extends the limited empirical basis 
of the MODIS aerosol algorithm io many environments found 
around the globe. 

The ratio of the reflectance n the red to that at 2.1 ym 
varies from slightly less than tht empirical value of 1/2 for 
dense and dark vegetation (p 2 ,i<0 1), to slightly more than 1/2 
for situations involving vegetation mixed with bright soils. This 
is in agreement with aircraft measurements over uniform dense 
vegetation [Renter et ai, 2001, Gatebe et a/., 2001] and over 
exposed soils but not for cases where vegetation is mixed with 
dark soils. One possible explanation for the differences for 
mixed scenes is the effect of shadows in an uneven terrain or 
uneven height of vegetation. Neither effect is explicitly included 
in the present model calculations. For example, a reflectance at 
2.1 ym of p 2 i=0.1 corresponds to a surface reflectance of 
p^O.035 in Fig. 2 or 3. However using values from Fig. 2 we 
can construct a surface cover with u combination of 60% of dark 
surface (Pr^O.06 and p red =0.015» with 40% of bright surface 
(p 2 i=0.30 and p red =0.18). If this combination is 35% shadowed 
and if these shadows are totally black in both channels, the 
linearly composited reflectance gri es the same reflectance at 2.1 
(p 2 !=0.10) but higher reflectance n the red of pr^pO.052, than 
the original value of p red =0.035. 

Nevertheless, and in agreement with the mentioned aircraft 
measurements, the conclusion of this study is that for very dark 
surfaces (e g., p 2 .i ^ 0.06) the surface reflectance in the red 
channel, may be overestimated by assuming a ratio of 1/2, 
resulting in an underestimate in the derived aerosol optical 
thickness. A smaller ratio shoi Id be used however more 
extensive field and modeling experiments have to be conducted 
in order to better quantify this ratic 
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